INITIAL DATA 

Material properties: porosity, water adsorption isotherm, 
pore solution, solid phase content, 
water content, diffusion coefficients, 
water diffusivity, ion adsorption 
isotherm, temperature 

Numerical parameters: number of lime steps 

length of time steps, FEM mesh 

Boudary conditions: concentrations, water content, 
electrical potential 
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TRANSPORT STEP 

Resolution of equations (2), (7), (14) with the finite 
element method 

OUTPUT: concentration profile for each ion 
water content profile electrical potential profile 
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CHEMICAL REACTION STEP 

Resolution of the chemical equilibrium algorithm 
at each node of the finite element mesh 
OUTPUT: corrected concentration for each ion 
and solid phase profiles 
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CHEMICAL DAMAGE 

Correction of the transport properties 
according to equation (44) 
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OUTPUT 

Solid phase profiles 
Ionic concentration profiles 
Adsorbed ion profiles 
later content profile 
Electrical potential profile 
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Applied electrical 
potential 



Agitating device 

40 

4- 



Upstream cell 

300 mmol/L NaOH 
500 mmol/L NaCI 



Downstream cell 
300 mmol/L NaOH 
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Saturated sample 
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Environment: 
70% Relative Humidity 
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Unsaturated soil: 85% relative humidity 
50 mmol/L S04 
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Electrical coupling and chemical activity — 
Electrical coupling without chemical activity - 
No electrical coupling nor chemical activity - 
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Paste sample 
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55 50 mmol/L of Na2S04 
in solution 
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Silicon coating 
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DETERMINING A FIRST CONCENTRATION 
OF IONS, A WATER CONTENT PROFILE AND 
AN ELECTRICAL POTENTIAL PROFILE USING 
A TRANSPORT ALGORITHM 



CALCULATING A CORRECTED CONCENTRATION ' 
OF IONS AND SOLID PHASE PROFILES USING A 
CHEMICAL REACTIONS ALGORITHM 




CACULATING A CHANGED 
TRANSPORT PROPERTIES PROFILE 



DETERMINING AN ION CONCENTRATION 
AND SOLID PHASE PROFILES 



INITIAL DATA 

Material properties: porosity, 

pore solution, 
temperature 
Numerical parameters: number of time steps 

length of time steps, FEM mesh, 
set of tortuosity 
Boudary conditions: concentrations 

electrical potential 



TRANSPORT STEP 

Resolution of equations (2), (14) with the finite 
element method with the diffusion coefficients 
corresponding to the tortuosity 
OUTPUT: concentration profile for each ion 
electrical potential profile 
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CALCULATE CURRENT 

Calculate the electrical 
current with equation (55) 
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OUTPUT 

The numerical current that 
match the measured current 
gives the diffusion 
coefficients for this 
material. 



DETERMINING AN ION CONCENTRATION 
AND AN ELECTRICAL POTENTIAL PROFILE 
USING A TRANSPORT ALGORITHM 
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DETERMINING A DIFFUSION COEFFICIENT 
FOR EACH ION USING THE CONCENTRATION AND 
THE ELECTRICAL POTENTIAL PROFILE 
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CALCULATING AN ELECTRICAL CURRENT 
USING TORTUOSITY PARAMETERS, , 
THE ION CONCENTRATION AND 
THE ELECTRICAL POTENTIAL 



82 



CHOOSING AN ELECTRICAL CURRENT 
CLOSEST TO A MEASURED ELECTRICAL CURRENT J 
AND DETERMINING A TORTUOSITY PARAMETER 
CORRESPONDING TO THE CHOSEN 
ELECTRICAL CURRENT 
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DETERMINING THE DIFFUSION COEFFICIENT 
OR EACH ION USING THE A 
CHOSEN TORTUOSITY 
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